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ABSTRACT 

We apply the global mode analysis, which has been recently developed for the modeling of kHz quasi-periodic 
oscillations (QPOs) from neutron stars, to the inner region of an accretion disk around a rotating black hole. 
Within a pseudo-Newtonian approach that keeps the ratio of the radial epicyclic frequency k to the orbital 
frequency fl the same as the corresponding ratio for a Kerr black hole we determine the innermost disk region 
where the hydrodynamic modes grow in amplitude. We find that the radiation flux emerging from the inner disk 
has the highest values within the same region. Using the flux weighted averages of the frequency bands over 
this region we identify the growing modes with highest frequency branches + k and to be the plausible 
candidates for the high-frequency QPO pairs observed in black hole systems. The observed frequency ratio 
around 1 .5 can therefore be understood naturally in terms of the global free oscillations in the innermost region 
of a viscous accretion disk around a black hole without invoking a particular resonance to produce black hole 
QPOs. Although the frequency ratio {n + K)/{n) is found to be not sensitive to the black hole's spin which 
is good for explaining the high-frequency QPOs it may work as a limited diagnostic of the spin parameter to 
distinguish black holes with very large spin from the slowly rotating ones. Within our model we estimate the 
frequency ratio of a high-frequency QPO pair to be greater than 1 .5 if the black hole is a slow rotator. For fast 
rotating black holes, we expect the same ratio to be less than 1 .5. 

Keywords: accretion, accretion disks — black hole physics — stars: oscillations — X-rays: stars 



1. INTRODUCTION 

The high energy emission from neutron stars and black 
holes in X-ray binaries is generally powered by accretion onto 
the compact object. The variability of X-ray light curve with 
different time scales from milliseconds to days is usually at- 
tributed to various characteristic time scales associated with 
accretion flow around the black hole or neutron star In low- 
mass X-ray binaries (LMXBs), where the central object is fed 
by an accretion disk, a few variability frequencies observed 
as quasi-periodic oscillation (QPO) peaks in addition to other 
broad-band features in the power spectra are common to both 
black hole and neutron star sources. Although there are some 
phenomenological differences between QPOs in black hole 
candidates and those observed in neutron star LMXBs, the 
similarities such as tight correlations of high and low fre- 
quency power spectral features in black hole and neutron star 
sources are remarkable (van der Klis 1994; Psaltis, Belloni, & 
van der Klis 1999; Wijnands & van der Klis 1999). Any inter- 
pretation or model solely based on the existence of a magnetic 
field, a hard surface, which are neutron star like properties, or 
an innermost stable circular orbit (ISCO) as a black hole like 
property to produce QPOs cannot account for the correlations 
of timing properties among different sources. 

QPOs were discovered in the X-ray power density spec- 
trum of black hole transients with frequencies in the ~ 0. 1 -30 
Hz range during the very high spectral state of these sources 
(Motch et al. 1983; Miyamoto & Kitamoto 1989; Miyamoto 
et al. 1991; Morgan, Remiflai-d, & Greiner 1997; Wij- 
nands, Homan, & van der Klis 1999; Sobczak et al. 2000; 
Strohmayer 2001; Muno et al. 2001). For several black hole 
sources, high-frequency QPOs were observed as single peaks 
roughly around 200 Hz (see Remillard & McClintock 2006 
and references therein). The discovery of twin kHz QPOs 
in neutron star LMXBs with their peak frequencies within 
around the 200- 1200 Hz range and a peak separation in the 

m.erkut@iku.edu.tr 



~ 180-360 Hz range together with their tight correlations 
with the low frequency power spectral components which 
are also observed in black hole candidates strengthens the 
idea that QPOs are produced in the inner regions of accre- 
tion disks around compact objects (see van der Klis 2000 and 
references therein; Mendez & Belloni 2007). The discover- 
ies of twin hectoHz QPOs from black hole candidates such 
as GRO J 1655-40, XTE J 1550-564, and GRS 1915+105 
(Remillard et al. 2002, 2003; Remillard & McCUntock 2006) 
have almost certified the idea of unifying the interpretation of 
high-frequency QPOs observed in black hole and neutron star 
sources within a single QPO model. 

One of the most striking differences between the high- 
frequency QPOs of black holes and the kHz QPOs in neutron 
star LMXBs is the fact that the former do not show any signif- 
icant correlation with the X-ray luminosity whereas the latter 
do. From one observation to another in a given source, the 
frequency shifts in the high-frequency QPOs of black holes 
are negligible as compared to the variations in the frequen- 
cies of kHz QPO peaks of accreting neutron stars. As com- 
pared to kHz QPOs detected in the power spectra of neutron 
stars, high-frequency QPOs from black holes are weak fea- 
tures having relatively low-quality factors. The commonly 
observed property of both low and high-frequency QPOs in 
black holes is the fact that these oscillations are strongest at 
photon energies above 6 keV when the power law compo- 
nent in the energy spectra dominates over the disk component 
(Remillard & McClintock 2006). Though it is not conclu- 
sive regarding the number of black hole sources exhibiting 
twin QPOs, the frequency ratio of the upper high-frequency 
QPO to the lower one is close to 3:2 in black hole candidates. 
There is no specific value for the ratio of two simultaneous 
kHz QPOs observed in neutron star sources; the ratio of the 
upper QPO frequency to the lower one rather takes different 
values changing between one and three from one source to 
another (Belloni, Mendez, & Homan 2005). Beside similari- 
ties and tight correlations, the phenomenological differences 
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between the high-frequency QPOs observed from black hole 
candidates and those from neutron star LMXBs likely arise 
from the dominant effect of different boundary conditions im- 
posed by a hard surface or a magnetic field and the ISCO as 
they might be more appropriate for a neutron star and a black 
hole, respectively. 

High-frequency QPOs were detected in seven black hole 
sources among which three black hole candidates exhibited 
QPO pairs in their power spectra. In all the sources with QPO 
pairs of commensurate frequencies, the ratio of two QPO fre- 
quencies is very close to 1 .5 in two black hole binaries and to 
1.6 in the third one (Remillard et al. 2002; Remillard et al. 
2003). The common property of high-frequency QPOs ob- 
served in black hole systems is based on the spectral state of 
a given source. All high-frequency black-hole QPOs are usu- 
ally observed in very high and high spectral states which are 
characterized by the most luminous states of the source (see 
Remillard & McClintock 2006 and references therein). For 
such high luminosities, the accretion disk around the black 
hole is expected to be truncated at the radius of the ISCO 
according to the standard model (Novikov & Thorne 1973; 
Shakura & Sunyaev 1973). For radii less than the radius of the 
ISCO, the accreting gas is thought to plunge radially towards 
the black hole. In the recent MHD simulations by Beckwith, 
Hawley, & Krolik (2008), however, the innermost ring of the 
disk which emits significant radiation has been shown to lie 
inside where the standard model predicts. These simulations 
have modified the stress-free boundary condition of the stan- 
dard model at the ISCO, but they have employed the classical 
relationship between magnetic stresses and energy dissipation 
and the assumption that the radial inflow time-scale of the ac- 
creting matter inside the ISCO is longer than the time-scales 
for thermalization and radiation of the dissipated heat, both 
of which may not be valid at all within the plunging region. 
We therefore anticipate the innermost part of the disk beyond 
the ISCO to be the most probable region near the black hole 
for the production of high-frequency QPOs. It is very likely 
within this region that the frequencies of the unstable growing 
disk modes correspond to the frequency bands of QPOs. 

The early attempts to interpret QPO frequencies in terms of 
disk modes were made by Alpar et al. (1992) and Alpar & Yil- 
maz (1997). The initial contributions to the theory of disk os- 
cillations explored the role of trapped disk oscillations in the 
variability of the X-ray power spectra of black hole sources 
(Kato & Fukue 1980; Nowak & Wagoner 1991; Kato 2001 
and references therein). The general relativistic test-particle 
frequencies were recognized and employed to construct mod- 
els of QPOs for both neutron stars and black holes (Stella, Vi- 
etri, & Morsink 1999; Abramowicz et al. 2003). The model 
of Psaltis & Norman (2000) revealed the importance of hy- 
drodynamic corrections to relativistic test-particle frequencies 
and the effect of hydrodynamic disk parameters on the corre- 
lations of QPOs and broad-band noise component. Alpar & 
Psaltis (2008) noted that the radial epicyclic frequency would 
be the highest dynamical frequency in the inner region of an 
accretion disk. This conclusion is based upon the existence of 
a magnetohydrodynamic boundary region around the neutron 
star where orbital frequencies deviate from Keplerian test- 
particle frequencies due to viscous and magnetic stresses (see 
Erkut & Alpar 2004). The recent analysis by Erkut, Psaltis, & 
Alpar (2008) of global hydrodynamic modes in the boundary 
regions of neutron stars showed how important the hydrody- 
namic effects are in estimating the observational characteris- 
tics of high-frequency QPOs. 



In this paper we apply the mode analysis, which has been 
developed by Erkut, Psaltis, & Alpar (2008) for the boundary 
region model of kHz QPOs from neutron stars, to the inner 
region of an accretion disk around a Kerr black hole. In order 
to account for the general relativistic effects of a Kerr metric 
on the ratio of dynamical frequencies, we work with a new 
pseudo-Newtonian potential that is appropriate for the analy- 
sis of temporal behavior of relativistic disks. This approach 
allows us to extend our recent study of global hydrodynamic 
modes of free oscillations to black holes as well. Our aim is to 
identify the modes whose frequency bands correspond to the 
high-frequency QPO pairs usually detected with a frequency 
ratio close to 1.5 in black hole sources. Most importantly, we 
provide a way to use the frequency ratio of these modes as a 
diagnostics of the spin parameter of a rotating black hole that 
exhibit high-frequency QPOs. 

In Section 2 we introduce our pseudo-Newtonian approach. 
The basic equations and parameters related to the analysis of 
disk modes for a Kerr black hole are presented in Section 3. 
In Section 4 we come up with the mode analysis and identify 
the modes relevant to the commensurate high-frequency QPO 
pairs with a frequency ratio around 1 .5. We discuss the results 
and present our conclusions in Section 5. 

2. PSEUDO-NEWTONIAN TREATMENT OF FREQUENCIES 

Using pseudo-Newtonian potentials in hydrodynamic sim- 
ulations is an effective and easy method to incorporate rel- 
ativistic effects into accretion flows (e.g., Chan, Psaltis, & 
Ozel 2009). The pseudo-Newtonian potential proposed by 
Paczyriski & Wiita (1980) is successful in estimating the 
ISCO for Schwarzschild geometry. It is therefore appropriate 
for a non-rotating black hole. Recently, two modified Newto- 
nian force models have been introduced by Mukhopadhyay & 
Misra (2003) to approximate the dynamical frequencies of an 
accretion disk around a rotating black hole. 

In our current analysis, the ratio of two successive fre- 
quency bands of disk modes is of interest in identifying the 
high-frequency QPOs from black holes. Unlike the early 
studies we mentioned above our pseudo-Newtonian treatment 
of dynamical disk frequencies keeps the ratio of the radial 
epicyclic frequency k to the orbital frequency 57 exactly the 
same as the corresponding ratio observed by a distant observer 
of a Kerr black hole. 

The Kerr expression for the ratio of the test-particle fre- 
quencies K and n is 

where a is the spin parameter of the black hole and R„ = 
GM/c^ with M and c being the mass of the black hole and 
the speed of light, respectively. The Newtonian expression 
for the same ratio is 

K r~f d\nVt\ 

The ratio given by equation (|2]i is valid for both the test- 
particles and the hydrodynamical fluids rotating in orbits. In 
the steady state of a geometrically thin disk, the hydrodynam- 
ical effects of pressure gradients on the test-particle frequen- 
cies are negligible. As we mention in Section 3, the radial 
momentum balance in a geometrically thin disk can be ap- 
proximated by the test-particle orbits where the centripetal ac- 
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celeration of each gas particle rotating with the frequency Q 
is due to the gravitational force. Setting equations ([T|) and (|2|i 
equal to each other, we obtain within our pseudo-Newtonian 
approach a differential equation for the orbital frequency fl{r). 
Given a suitable pseudo-Newtonian potential, the orbital fre- 
quency fl satisfies the radial momentum equation and mimics 
the effect of strong gravity on the gas-particle orbits by keep- 
ing the value of n/fl the same as the corresponding Kerr value 
in the test-particle regime. In this sense, our approach is sim- 
ilar to the early pseudo-Newtonian treatments of the dynami- 
cal frequencies in a geometrically thin accretion disk around 
a black hole. To find a solution for ri(r), we require that our 
pseudo-Newtonian orbital frequency match the Kerr orbital 
frequency in the outer disk. According to a distant observer, 
the Kerr expression for the orbital frequency is 



1+fl {Rg/r 



3/2 ' 



(3) 



1/2 



where il.K{r) = [GM / r^y " is the Keplerian frequency. For 
sufficiently large radii, that is for r ^ R„, the Kerr orbital fre- 
quency assumes its Keplerian value. Using the same asymp- 
totic boundary condition on the pseudo-Newtonian orbital fre- 
quency, it follows from equations ([T]i and (|2]i that 



n{r) = i7K('")exp 



8 /R 



3/2 



+ -fl 

4 



R. 



(4) 

To illustrate our treatment of the radial epicyclic and or- 
bital frequencies, we plot in Figure 1 the pseudo-Newtonian 
frequencies k and fl over a wide range of disk radii in com- 
parison with the corresponding Kerr frequencies. Figure la 
shows the radial profiles of k and 51 for a Schwarzschild black 
hole for which the spin parameter = 0. The radial profiles of 
the same frequencies for a Kerr black hole with a spin param- 
eter a = 0.8 are shown in Figure lb. The pseudo-Newtonian 
frequencies k and ft can be seen to slightly deviate from the 
Kerr frequencies in the inner disk while they asymptotically 
match them in the outer disk. Note, however, that the pseudo- 
Newtonian frequency k always yields the correct estimation 
for the ISCO as n/ft matches its Kerr value exactly for all 
disk radii. 

In the next section, we present the basic equations for a 
geometrically thin disk and obtain within the current pseudo- 
Newtonian approach the hydrodynamic parameters that are 
necessary for the analysis of global modes in the inner disk. 

3. BASIC EQUATIONS AND PARAMETERS 

The long-wavelength global hydrodynamic modes of free 
oscillations have been recently studied for the boundary re- 
gions of accretion disks around neutron stars (see Erkut, 
Psaltis, & Alpai- 2008, hereafter EPA08). In the mode analy- 
sis by EPA08 the basic disk equations are perturbed for a geo- 
metrically thin disk in vertical hydrostatic equilibrium. The 
Fourier decomposition of the linearized perturbation equa- 
tions leads to the identification of the complex mode frequen- 
cies The real parts of axisymmetric (m = 0) and nonax- 
isymmetric (m > 1) mode frequencies correspond to the fre- 
quency bands of QPOs while the imaginary parts determine 
the growth rates of the oscillations. 

In the global mode analysis, both the frequency bands and 
their growth rates depend on several key parameters such as 
n/ft, the radial profile of the surface density, (3, the inverse 



timescale fl^, associated with the radial drift velocity, and the 
inverse timescale fJ, associated with the sound speed in the 
inner disk. These parameters are determined by the global 
structure of the unperturbed steady disk (see EPA08). 

For sufficiently high mass accretion rates, e.g., M > 0. IMe, 
where Me is the Eddington mass accretion rate, black holes 
in LMXBs accrete matter through radiatively efficient accre- 
tion disks whose innermost regions are dominated by radia- 
tion pressure (Shakura & Sunyaev 1973, hereafter SS73). The 
innermost truncation radius of such a disk around a Kerr black 
hole is estimated by the radius of the ISCO, rin, which can be 
found as a solution of K/fi = for r = rin (see eq. IT]). The un- 
perturbed steady structure of a radiation pressure dominated 
inner disk is described by 



PC. = 3 , 



(5) 



where is the effective sound speed, p is the average mass 
density, and e is the radiation energy density. We write, for 
the vertical hydrostatic equilibrium in the disk. 



(6) 



where H is the half-thickness of the disk. The average mass 
density can be expressed in terms of the surface mass density 
S as 

E 



(7) 



The vertical energy balance in the inner disk is satisfied for 



4c 



(8) 



where Kes is the electron scattering opacity and $ is the en- 
ergy dissipation rate per unit area of the disk (see SS73). The 
energy flux due to viscous energy dissipation is 



$ = -i^S r 

2 \ dr 



(9) 



Here, i/ is the kinematic viscosity for which the a-prescription 
(SS73) can be written as 



dr 



(10) 



For a geometrically thin disk, the radial momentum balance 
can be written to a good approximation as fl-r-dT/dr = 0, 
where F is the pseudo-Newtonian potential that mimics the 
gravitational field of the black hole. In the radial momentum 
equation, the pseudo-Newtonian force, -dT /dr, is the source 
of acceleration, -il^r, where O is given by equation (|4]i. For 
the conservation of mass and angular momentum, we write 



and 



-2T:rY,Vr=M 



-, dil ■ o 

27rzyl]r^ +Mr^n = C, 

dr 



(11) 



(12) 



respectively, where v,- is the radial drift velocity of the accret- 
ing matter in the inner disk and C is an arbitrary constant of 
integration. 

We solve equation ST% using torque-free boundary condi- 
tion at the innermost disk radius, ^n, which is appropriate for 
a disk around a black hole. The constant of integration can be 
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determined as C =Mrf^^il{ri„) to satisfy the torque-free bound- 
ary condition. Using equation (|4|i, it follows from equations 
©-([nil that 



flflnS 3 ^ Dir) E(r) 
lilnr 2 A(r) F(r) 



EE ~ 2.57am — ^ A(r), 

Cs/r \ r J 

-^-2.4mf — ^jA(r)F(r), 



and 



f{r)A(r)F(r), 



(13) 
(14) 
(15) 

(16) 



where $t = 3GMM/87rr?j is the typical value for the radiation 
flux and 



M 



M (jL)(M.y (17) 

IQi^gs-i Vo.06/ \MqJ 



(18) 



Me 1.9xlOi«g5 
Here, 77 is a function of the spin parameter a such that 

^^_ 4V2[l-(l-r;)^]'^'-2(l-77) 

3V3[l-(l-77)'] 

for prograde accretion disks around rotating black holes (see 
Shapiro & Teukolsky 1983). In equations the di- 

mensionless factors arising from the boundary conditions and 
pseudo-Newtonian corrections are 



B(r) = -3 ( ^ l+4fl( 



3/2 



Z)(r) = -4('^')+8fl (^^''' 



3 2(R, 



4a' f ^ 



-1/2 



and 



F(r)=l- 



fir) \n 



-1/2 



with 

f(r) = exp 



R„\ 8 



3/2 



(19) 
(20) 
(21) 
(22) 

(23) 
(24) 



For illustrative purposes we display in Figure 2 the radial 
distributions of the outgoing radiation flux (see eq. |[T6l ) 
throughout the inner disk for two putative black holes with 
spin parameters a = and a = 0.9. In the following section, 
we use the ratio n/il (see eq. HI) and the global hydrody- 
namic parameters /?, il^/fls, and njn (see eq. iflJll - fTSl ) 
to identify the radial zone in the inner disk where the modes 
grow. As we will see, the hydrodynamic modes grow only 
within a limited range of radii in the innermost disk region 
out of which the radiation flux is maximum (see Fig. 2). 



4. GLOBAL MODES IN THE INNER DISK 

When there are no external perturbations due to large-scale 
magnetic fields of the accreting star, the free oscillation modes 
in a boundary region or the inner disk are excited through the 
dynamical effect of the viscosity. In the limit of small hydro- 
dynamic corrections, this can be seen from the growth rates 
of both axisymmetric and nonaxisymmetric modes for which 
Im(tj'"'*) cx -prii, (see EPA08). In the presence of viscos- 
ity and therefore of radial drift velocity, 57^ ^ and the high- 
frequency modes can have positive growth rates only if /3 < 0. 
This is also valid for the global modes in a disk around a black 
hole. Unlike neutron stars, the effect of a large-scale toroidal 
magnetic force in addition to viscosity on the excitation of 
global modes might be absent in the case of black holes (see 
Section 3 in EPA08). 

For a black hole disk, the presence of the ISCO with a 
torque-free boundary condition determines the radial profile 
of the surface density, /3, and thus the growth rates of the 
modes. Note that ^ ~ 3/2 > for r > (see eq. |fT3l ) 
and we expect, in the regime of small hydrodynamic correc- 
tions, that the global hydrodynamic modes do not grow for 
sufficiently large radii in the inner disk. For the innermost 
disk region, however, the hydrodynamic corrections can be 
important to distinguish among the growth rates of different 
modes and to identify the set of radii at which these modes 
grow. In order to see the effects of hydrodynamic parameters 
on both the frequency bands and growth rates of the modes in 
the inner disk beyond the regime of negligible hydrodynamic 
corrections, we use equations (fT3Tl-(fT5Tl together with equa- 
tion ([T) in the full eigenfrequency solutions for axisymmetric 
and nonaxisymmetric perturbations given in the Appendix of 
EPA08. 

Figures 3-6 show the real and imaginary parts of the com- 
plex mode frequencies in units of the orbital frequency 57 as 
functions of the radial distance in the inner disk. The real 
and imaginary parts represent the frequency bands and the 
growth rates of the modes, respectively. In Figures 3-6, we 
label the hydrodynamic mode frequencies and their growth 
rates with notation corresponding to the test-particle frequen- 
cies. This provides us with an easy identification and a sim- 
ple designation of each mode without ambiguity and without 
loss of generality. In the limit of small hydrodynamic correc- 
tions, the frequencies of all hydrodynamic modes converge to 
the test-particle frequencies. We mark axisymmetric (m = 0) 
modes with the corresponding test-particle frequencies uj = 
and uj = K. We use uj = 0, and w = ± k as the appropriate 
labels to distinguish among nonaxisymmetric modes {m = 1). 

Figure 3a exhibits the run of the mode frequencies in the 
inner disk of a Schwarzschild black hole (a = 0). We display 
the growth rates of the modes in Figure 3b. Figure 3 is ob- 
tained for a = 0.1 and rii = Q).\. Note that the hydrodynamic 
modes do not grow in the inner disk for r > 2.3ri„. The ra- 
dial zone within which all modes grow covers only a limited 
range of radii around - 2 . 1 rin in the innermost disk region, 
as shown in Figure 3b. The hydrodynamic modes with fre- 
quency bands around il and have relatively higher growth 
rates as compared to those around il ± k and n bands. Fig- 
ure 4 reveals how the mode frequencies and the growth rates 
are affected by the rotation of the black hole. For the same 
values of the viscosity parameter and the mass accretion rate, 
that is, for a = 0. 1 and ;7i = 0. 1 , we plot the frequency profiles 
of the modes in Figure 4a and the corresponding growth rates 
in Figure Ab for a rotating black hole with a spin parameter 



ON THE HIGH-FREQUENCY QPOs FROM BLACK HOLES 



5 



a = 0.9. In comparison with Figure 3b, the growth rates of the 
modes are higher in Figure 4b. The range of radii at which the 
hydrodynamic modes grow is around ri„-2.3rin. All modes 
decay for r > 2.5ri„ (see Fig. 4b). We obtain Figure 5 and Fig- 
ure 6 keeping the spin parameter of the black hole at a = 0.9, 
however, changing the viscosity parameter a and the mass ac- 
cretion rate m. For a = 0.01 and m = 0. 1, we explore the run of 
the mode frequencies and the growth rates in Figures 5a and 
5b, respectively. As compared to Figure 4a, the radial profiles 
of the mode frequencies can be seen to be almost unaffected 
by a change in the viscosity parameter a (see Fig. 5fl). We ob- 
serve, in Figure 5b, that the growth rates are lower than those 
in Figure 4b by a factor around 0. 1 which, indeed, is the factor 
of decrease in a. Figure 6 is obtained for m = 0.6 while keep- 
ing the values of a and the spin parameter a the same as in 
Figure 5. Note that both the frequency bands and the growth 
rates of the modes are modified to some level at relatively high 
mass accretion rates. The greater the mass accretion rate m, 
the higher are the growth rates of hydrodynamic modes (see 
Fig. 6b). We note that the frequency bands that are related to 
fl±K, fi, and K branches in the limit of small hydrodynamic 
corrections begin to deviate from the test-particle frequencies 
for sufficiently large mass accretion rates (jn > 0.6) as shown 
in Figure 6a. 

The common property of Figures 3-6 is that all the hydro- 
dynamic modes grow within a limited region in the innermost 
part of the disk with characteristic radii in the ~ rin-2.5rin 
range. The radiation flux emerging from the same region 
attains the highest values with maxima at r ~ 1.6rin and 
r ~ 1 .Irin for the black holes with spin parameters a = and 
a = 0.9, respectively (see Fig. 2). It is interesting to deduce 
from Figures 3-6 that the frequency ratio of the hydrodynamic 
modes we associate with H+k and ft frequency bands is close 
to 1 .5 at radii in the ~ (L6- 1 .7)ri„ range, where the disk flux 
is maximum. As mentioned above, the modes, however, grow 
throughout an extended region, with < r < 2.5rin, of radii 
rather than being excited at a particular radius. Moreover, 
the modes with frequencies around 57 - k and k bands also 
grow within the same region. To distinguish among the pairs 
of growing modes which can be regarded as plausible candi- 
dates for the high-frequency QPO pairs from black holes, we 
consider the mutual ratios of the flux weighted averages of 
the frequency bands for different modes. We define the flux 
weighted average of a frequency branch uj as 

(uj) = [ 2Trruj(r)<S>(r)dr / f 2Tir^(r)dr , (25) 

J rm I J r,„ 

where is the critical radius beyond which the correspond- 
ing mode decays in the inner disk. Using equation ( l25T l. we 
calculate the ratios of the flux weighted averages of the fre- 
quency bands n + k, ft, ft- k, and k for different values of the 
black hole spin parameter a between and 1 . For each model 
value of a, we find a critical radius such that all the hy- 
drodynamic modes grow for r < r^. Unlike the growth rates, 
the mode frequencies and the width of the radial zone where 
the modes grow are sensitive to the spin parameter a, but not 
sensitive to the viscosity parameter a and the mass accretion 
rate m (see Figs. 3-6). In Figure 7 we display the run of 
+ {n)/{K), and {n)/{n-K.) forO <« < l. Figure? 

is obtained for the typical values, a = 0.1 and m = 0.1. The 
values for {ft + are densely clustered around 1 .5 over a 

wide range of values for a as shown in Figure 7. For slow ro- 
tators (a ~ 0), {ft + K)/{ft) ~ 1 .6. The same ratio drops below 



1.5 as fl approaches 1 for rapidly rotating black holes. The 
values of (57)/(k) and k), on the other hand, span a 

wide range as the spin parameter a varies between and 1 
(see Fig. 7). Our analysis suggests the hydrodynamic modes 
with frequency bands around n + k and 57 to be the plausi- 
ble candidates for the high-frequency QPO pairs observed in 
black hole systems. Note that our model estimation for the 
frequency ratio of high-frequency QPO pairs involves the two 
highest frequency modes with positive growth rates. 

We give examples for surface density perturbations of such 
global modes in Figures 8-11. We display the three di- 
mensional profile of surface density perturbation Ei in terms 
of background surface density So in the innermost region 
(Hn <f'< 3nn) of a disk around a rotating black hole with spin 
parameter a = 0.9 for the typical values, a = 0.1 and m = 0.1. 
The examples for nonaxisymmetric modes with frequencies 
ft (Fig. 8) and ft+n (Fig. 9) show the surface density per- 
turbations at the time t = 2Pin, where Pi„ = 27r/i7(rin) is the 
rotation period at the innermost disk radius. In Figures 8 and 
9, the spiral like shapes of different iso-level contours plotted 
on the xy-plane reveal the similar nonaxisymmetric nature of 
these modes. In the long run, such as for t > lOOPi,,, the sur- 
face density perturbations of both axisymmetric and nonax- 
isymmetric modes grow in amplitude only for r < 2rin within 
the same domain. We illustrate this typical behavior in Fig- 
ure 1 1 as compared to Figure 10 for the case of axisymmetric 
mode with frequency k. Note that the perturbations at the 
time f = 15Pin (see Fig. 10) are comparable in amplitude over 
the whole computational domain (nn < r < 3rin). The pertur- 
bations at the time t = 200Pin, however, have large amplitudes 
only for < < 2rin whereas their amplitudes become neg- 
ligible for r > 2rin (see Fig. 1 1). 

In the global three-dimensional magnetohydrodynamic 
simulations of black hole accretion disks, the innermost disk 
region near the ISCO was found to show QPOs with fre- 
quency around the maximum of epicyclic frequency (Machida 
& Matsumoto 2003). In one of the recent simulations of 
the three-dimensional magnetohydrodynamic accretion flows 
around Schwarzschild black holes (Kato 2004), the structure 
of the flow has been changed at radial distances within the 
3.8 rs < r < 6.3 rs range, where rs is the Schwarzschild ra- 
dius. Two pairs of QPOs have been observed to be excited in 
that region with frequencies around the Keplerian frequency 
and the sum of Keplerian and epicyclic frequencies in the 
power spectra of these simulations. Most importantly, the fre- 
quency ratio of these QPO features has been found to be near 
1 .5. These results are in close agreement with the result of our 
mode analysis in the present work. 

5. DISCUSSION AND CONCLUSIONS 

We have probed the stability of the global modes in the in- 
ner region of a standard accretion disk around a black hole. 
Our study is the application of the recently developed analy- 
sis of global hydrodynamic modes (see EPA08) to the iden- 
tification of the high-frequency QPO pairs observed in black 
hole sources. The presence of the ISCO allows for effects 
of strong gravity on both the dynamical frequencies and the 
global hydrodynamic parameters. Our pseudo-Newtonian ap- 
proach takes account of these effects to determine the fre- 
quency bands and the growth rates of the unstable modes in 
the inner disk. 

The disk is truncated at the radius of the ISCO, ^n. The 
growth rates of the modes are negative for sufficiently large 
distances from the ISCO. We find that the modes grow in am- 
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plitude only within a narrow zone in the innermost disk re- 
gion. For a non-rotating black hole (a = 0) the characteristic 
radii of the zone lie in the ~ rin-2.1rin range. The modes 
grow within the ~ Hn — 2.5rin range for a rotating black hole 
with spin a = 0.9. Among the growing modes the growth 
rates of the frequency branches around il and are higher as 
compared to those of the modes with frequency bands around 
n±K and K (see Figs. 3-6). Due to the effect of enhanced hy- 
drodynamic corrections on the growth rates, the modes grow 
faster in an accretion regime with relatively high rate and vis- 
cosity (see EPA08). 

The radiation flux due to viscous energy dissipation in the 
inner disk takes the highest values within the narrow region 
where the modes grow (see Fig. 2). We deduce from the ra- 
dial profiles of the mode frequencies that the frequency ratio 
of the modes around fl + n and fl bands is very close to 1 .5 at 
the radius where the disk radiation is maximum. This value 
was observed for the frequency ratio of the high-frequency 
QPO pairs in black hole sources (see Remillard & McClin- 
tock 2006). Instead of being excited at a particular radius in 
the disk, the hydrodynamic modes grow in a region of finite 
radial extension. To make an estimation for the expectation 
value of a frequency band and therefore for the frequency ra- 
tios of the relevant modes, we calculate the flux weighted av- 
erages of the frequency bands over the innermost disk region 
where the modes grow. Scanning the ratios of the expected 
mode frequencies for all possible values of the spin parame- 
ter a (see Fig. 7), we find that only the modes around H + k 
and il branches have a frequency ratio around 1.5. This ra- 
tio is slightly higher than 1 .5 if the black hole is a slow ro- 
tator (O < fl < 0.8). The same ratio falls below 1.5 for fast 
rotating black holes (0.9 < a < l). The frequency ratios of 
other modes significantly deviate from 1 .5 over a wide range 
of values for the spin parameter. Relying on the observed val- 
ues for the frequency ratio of the upper high-frequency QPO 
to the lower one, we conclude that the modes with frequency 
branches H+k and il are the most plausible candidates for the 
high-frequency QPOs from black holes. 

The observations of high-frequency QPOs can be used to 
determine the underlying mechanism that produces these os- 
cillations and to measure the spin parameter of the black hole 
(Remillard & McClintock 2006). Our analysis may provide 
a way to employ the observed frequency ratio of a high- 
frequency QPO pair in a given source to estimate the spin 
parameter a. In this sense. Figure 7 comes out as an efficient 
tool for reading the spin parameter a that corresponds to the 
value of (ri + K)/(ri) to be interpreted as the frequency ratio 
of a high-frequency QPO pair observed in the X-ray power 
spectra of the black hole source. 

There are several reasons for expecting to see the finger- 
print of global long wavelength modes in the form of high- 
frequency QPOs observed in the X-ray power spectra of black 
hole sources in LMXBs. As compared to the neutron stars, 
there is little chance for accretion flows around the black holes 
in LMXBs to be affected by the dynamical action of a mag- 
netic field of stellar origin. In the case of a black hole, instead 
of a direct feedback from the compact object, except grav- 
ity, the fluctuations in the mass transfer rate from the binary 
companion introduce perturbations with a broad band of fre- 
quencies including those of the inner disk. The disk modes 
which depend on global disk parameters become unstable in 
the innermost disk region and thus the disk oscillation fre- 
quencies are selectively amplified without any need for an 



external mechanism to force them to attain high amplitudes. 
Furthermore, the observable luminosity variation in the X-ray 
light curve of a source due to global free oscillation modes 
of long wavelength, that is, of sufficiently large lengthscale 
is expected to be least affected by the MHD turbulent eddies 
of short wavelength. According to our present analysis, the 
higher the mass accretion rate M and the greater the viscos- 
ity parameter a, the higher are the growth rates of the modes. 
We therefore expect to observe these modes particularly in the 
state of high mass accretion rate and high viscosity. In such 
a state, the turbulent disk may also interact with its corona 
(see Tagger & Varniere 2006). This would lead to the for- 
mation of high-frequency QPOs in a spectral state where the 
contribution from the power law component is important. In 
our present analysis we identify the relevant disk modes with- 
out deliberating the disk-corona interaction which we plan to 
consider in a future work. 
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Figure 1. Radial profiles of the orbital and epicyclic frequencies. The Schwarzschild orbital and radial epicyclic frequencies fisch and ftsch and the Kerr orbital 
and radial epicyclic frequencies ^Kerr and ftKerr are shown by the solid curves. The pseudo-Newtonian orbital frequency Q and the pseudo-Newtonian radial 
epicyclic frequency k are shown by the dashed curves. The panel (a) is obtained for a non-rotating black hole with spin parameter a = 0. The frequencies in panel 
(b) are plotted for a rotating black hole with spin parameter a = 0.8. 
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Figure 2. Radiation flux emerging from the inner disk. In both the panels (a) and (b), the solid curve represents the general relativistic (GR) estimation for the 
disk flux according to Page & Thome (1974). The dashed curve in both the panels (a) and (b) estimates the disk flux according to our pseudo-Newtonian (PN) 
analysis. The panel (a) shows the flux from the accretion disk aroimd a non-rotating black hole (a = 0). The disk flux for a rotating black hole with spin parameter 
a = 0.9 is displayed in panel (b). 
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Figure 3. Real and imagmary parts of the complex frequencies for axisymmetric (m = 0) and non-axisymmetric (m = 1) modes in the innermost region of the 
accretion disk around a non-rotating black hole (a = 0). The viscosity parameter and the accretion rate are a = 0.1 and m = 0.1, respectively. The mode frequencies 
in panel (a) and the growth rates of the modes in panel (b) are labeled with the corresponding test-particle frequency branches. 
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Figure 4. Real and imaginary parts of the complex frequencies for axisymmetric (m = 0) and non-axisymmetric (m = 1 ) modes in the innermost region of the 
accretion disk around a rotating black hole with spin parameter a = 0.9. The viscosity parameter and the accretion rate are a = 0. 1 and m = 0.1, respectively. The 
mode frequencies in panel (a) and the growth rates of the modes in panel (b) are labeled with the corresponding test-particle frequency branches. 
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Figure 5. Real and imaginary parts of the complex frequencies for axisymmetric (m = 0) and non-axisymmetric (m = 1) modes in the innermost region of the 
accretion disk around a rotating black hole with spin parameter a = 0.9. The viscosity parameter and the accretion rate are a = 0.01 and m = 0.1, respectively. 
The mode frequencies in panel (a) and the growth rates of the modes in panel (b) are labeled with the corresponding test-particle frequency branches. 
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Figure 6. Real and imaginary parts of the complex frequencies for axisymmetric (m = 0) and non-axisymmetric (m = 1) modes in the innermost region of the 
accretion disk around a rotating black hole with spin parameter a = 0.9. The viscosity parameter and the accretion rate are a = 0.01 and rh = 0.6, respectively. 
The mode frequencies in panel (a) and the growth rates of the modes in panel (b) are labeled with the corresponding test-particle frequency branches. 
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Figure 7. Ratios of the flux weighted averages of different frequency bands of the modes growing in the iimer disk for all possible values of the spin parameter 
a between and 1. The frequency ratios as functions of the spin parameter are obtained for the typical values, a = 0.1 and m = 0.1. 




Figure 8. Profile of surface density perturbation S i with respect to background surface density So for nonaxisymmetric mode with frequency f2 in the innermost 
region (rin < < 3ri„) of a disk aroimd a rotating black hole with spin parameter a = 0.9 for the typical values, a = 0.1 and rh = 0.1. The profile is obtained at the 
time t = IPin, where Pm = In/Clinn) is the rotation period at the imiermost disk radius. 




Figure 9. Profile of surface density perturbation S i with respect to background surface density Sq for nonaxisymmetric mode with frequency Q + k in the 
innermost region (rm < r < 3rin) of a disk around a rotating black hole with spin parameter a = 0.9 for the typical values, a = 0.1 and m = 0.1. The profile is 
obtained at the time t = IPjn, where Pin = 27r/r2(ri„) is the rotation period at the innermost disk radius. 
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Figure 10. Profile of surface density perturbation S i with respect to background surface density So for axisymmetric mode with frequency k in the innermost 
region (rin < < 3ri„) of a disk around a rotating black hole with spin parameter a = 0.9 for the typical values, a -0.1 and m = 0. 1 . The profile is obtained at the 
time t = 15Pm> where Pm = 27r/n(rin) is the rotation period at the innermost disk radius. 




Figure 11. Profile of surface density perturbation Si with respect to background surface density Sq for axisymmetric mode with frequency k in the innermost 
region (n„ < r < Snn) of a disk around a rotating black hole with spin parameter a = 0.9 for the typical values, a = 0.1 and m = 0.1. The profile is obtained at the 
time t = 200Pin, where = 27r/f2(ri„) is the rotation period at the innermost disk radius. 



